function F=myfun1(x)

alpha=0.5;
sigma=5;
T=1.5;
mu=0.4;
global lambda1;

g1=x(1);
g2=x(2);
lambda2=x(3);
w2=x(4);
E1=x(5);
E2=x(6);


F=[g1^(1-sigma)-lambda1*g1^(-alpha*sigma)-lambda2*w2^(1-sigma*(1-alpha))*g2^(-sigma*alpha)*T^(1-sigma);
g2^(1-sigma)-lambda1*g1^(-alpha*sigma)*T^(1-sigma)-lambda2*w2^(1-sigma*(1-alpha))*g2^(-sigma*alpha);
g1^(alpha*sigma)/(1-alpha)-E1*g1^(sigma-1)-E2*g2^(sigma-1)*T^(1-sigma);
w2^(sigma*(1-alpha))*g2^(alpha*sigma)/(1-alpha)-E1*g1^(sigma-1)*T^(1-sigma)-E2*g2^(sigma-1);
E1-mu-alpha*lambda1/(1-alpha);
E2-mu*(w2*lambda2-lambda2+1)-alpha*w2*lambda2/(1-alpha)];